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Abstract 

We present the first measurement of the D* + width using 9/fb of e + e _ data 
collected near the T(4S) resonance by the CLEO II. V detector. Our method 
uses advanced tracking techniques and a reconstruction method that takes 
advantage of the small vertical size of the CESR beam spot to measure the 
energy release distribution from the D* + — > D -k + decay. We find T(D* + ) = 
96 ± 4 (Statistical) ± 22 (Systematic) keV. We also measure the energy release 
in the decay and compute Am = m D *+ — m^o = 145.412±0.002(Statistical) ± 
0.012 (Systematic) MeV/c 2 . 
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A measurement of T(D* + ) opens an important window on the non-perturbative strong 
physics involving heavy quarks. The basic framework of the theory is well understood, 
however, there is still much speculation - predictions for the width range from 15keV to 
150 keV 0. We know the D* + width is dominated by strong decays. The level splitting in 
the B sector is not large enough to allow real strong transitions. Therefore, a measurement 
of the width of the D* + gives unique information about the strong coupling constant in 
heavy-light meson systems. This width only depends on g, a universal strong coupling 
between heavy vector and pseudoscaler mesons to the pion, since the small contribution of 
the electromagnetic decay can be neglected, yeilding 

where f n is the pion decay constant and the momenta are for the indicated particle in D* + 
decay in the D* + rest frame 0. 

Prior to this measurement, the D* + width was limited to be less than 131 keV at the 90% 
confidence level by the ACCMOR collaboration |5). This letter describes a measurement of 
the D* + width with the CLEO II. V detector ||J]. The signal is reconstructed through a single, 
well-measured sequence, D* + — > ty^ ow D°, D° — > K~n + . Consideration of charge conjugated 
modes are implied throughout this letter. 

The CLEO detector has been described in detail elsewhere. All of the data used in this 
analysis are taken with the detector in its II. V configuration ||. The data were taken in 
symmetric e + e~ collisions at a center of mass energy around 10 GeV with an integrated 
luminosity of 9.0/fb provided by the Cornell Electron-positron Storage Ring (CESR). The 
nominal sample follows the selection of D* + — > 7i^ ow D° — > K~ir + 7r^ low candidates used in our 
D° — D° mixing analysis ||. 

Our reconstruction method takes advantage of the small CESR beam spot and the kine- 
matics and topology of the D* + — > tt^ ow D° — > 7r^ ow K~n + decay chain. The K~ and tc + are 
required to form a common vertex. The resultant D° candidate momentum vector is then 
projected back to the CESR luminous region to determine the D° production point. The 
CESR luminous region has a Gaussian width ~ 10 /zm vertically and ~ 300 /im horizontally. 
This procedure determines an accurate D° production point for D 0, s moving out of the hor- 
izontal plane; D 0, s moving within 0.3 radians of the horizontal plane are not considered. 
Then the 7r^ ow track is refit constraining its trajectory to intersect the D° production point. 
This improves the resolution on the energy release, Q = M{K~^ + ^ Xow ) — M(K~tx + ) —m w +, 
by more than 30% over simply forming the appropriate invariant masses of the tracks. The 
improvement to resolution is essential to our measurement of the width of the D* + . The 
distribution of our resolution, erg, is shown in Figure [l] and is typically 150 keV. The good 
agreement between Monte Carlo and data demonstrates that the kinematics and sources 
of uncertainties on the tracks, such as the number of hits used and the effects of multiple 
scattering in detector material, are well modeled. 

The challenge of measuring the width of the D* + is understanding the tracking system 
response function since the experimental resolution exceeds the width we are trying to mea- 
sure. We depend on exhaustive comparisons between a GEANT based detector simulation 
and our data. We addressed the problem by selecting samples of candidate D* + decays using 
three strategies. 
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FIG. 1. Distribution of ctq, the uncertainty on Q as determined from propagating track fitting 
errors. The arrow indicates a selection discussed in the text. 
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First we produced the largest sample from data and simulation by imposing only basic 
tracking consistency requirements. We call this the nominal sample. Second we refine the 
nominal sample selecting candidates with the best measured tracks by making very tight 
cuts on tracking parameters. We call this the tracking selected sample. A third alternative 
is to select our data on specific kinematic properties of the D* + decay that minimize the 
dependence of the width of the D* + on detector mismeasurements. We call this the kinematic 
selected sample. In all three samples the width is extracted with an unbinned maximum 
likelihood fit to the energy release distribution and compared with the simulation's generated 
value to determine a bias which is then applied to the data. These three different approaches 
yield consistent values for the width of the D* + giving us confidence that our simulation 
accurately models our data. 

To further improve the quality of reconstruction in our sample, we apply some selections 
at the kinematic boundaries of 7r^ ow momentum and the opening angle 8 between the vr^" ow 
and the D° candidate as a function of the D* + candidate momentum distributions to remove 
a small amount of misreconstructed signal and background. We also require o~q < 200 keV 
which removes the long tail in the error distribution. 

Table [D] summarizes the statistics in our three samples. The tracking and kinematic 
samples are subsets of the nominal sample. The two subsets contain 94 common candidates. 

We assume that the intrinsic width of the D° is negligible, T(D°) <C T(D* + ), implying 
that the width of Q is simply a convolution of the shape given by the D* + width and 
the tracking system response function. Thus we consider the pairs of Q and oq for D* + — > 
7r^" ow D° — > K~7T + 7r^ ow where o~q is given for each candidate by propagating the tracking errors 
in the kinematic fit of the charged tracks. We perform an unbinned maximum likelihood fit 
to the Q distribution. 

The underlying signal shape of the Q distribution is assumed to be given by a P-wave 
Breit-Wigner with central value Qo- We considered a relativistic and non-relativistic Breit- 
Wigner as a model of the underlying signal shape, and found negligible changes in the fit 
parameters between the two. The width of the signal Breit-Wigner depends on Q and is 
given by 



where r = T(D* + ), P and M are the candidate vr^ ow or D° momentum in the D* + rest 
frame and iT7T7r slow mass, and P and M are the values computed using Q . The effect of the 
mass term is negligible at our energy. The partial width and the total width differ negligibly 
in their dependence on Q for Q > 1 MeV. 

For each candidate the signal shape is convolved with a resolution Gaussian with width 
<tq, determined by the tracking errors, as a model of our finite resolution shown in Figure |I[ 
The fit also includes a background contribution with a fixed shape derived from our simu- 
lation, and modeled with a third order polynomial. We allow a small fraction of the signal, 
f m i s , to be parametrized by a single Gaussian resolution function of width a m i S . This shape 
is included in the fit to model the tracking mishaps which our simulation predicts to be at the 
5% level in the nominal sample and negligible in both the tracking and kinematic selected 
samples. In our standard fit we constrain the level of this contribution while allowing o- mis 




(2) 



to float. 
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FIG. 2. Fits to the three data samples: a) is nominal; b) is tracking; and c) is kinematic. The 
different contributions to the fits are shown by different colors. 



TABLE I. Results of the fits described in the text. The uncertainties are statistical. 





Sample 


Parameter 


Nominal 


Tracking 


Kinematic 


To (keV) 


98.9 ± 4.0 


106.0 ± 19.6 


108.1 ±5.9 


Qo (keV) 


5853 ± 2 


5854 ± 10 


5850 ± 4 


N s 


11207 ± 109 


353 ± 20 


3151 ±57 


fmis (%) 


5.3 ±0.5 


NA 


NA 


a mis (keV) 


508 ± 39 


NA 


NA 


N b 


289 ± 31 


15 ±7 


133 ± 16 



As a preliminary test to fitting the data we run the complete analysis on a fully simu- 
lated sample that has about ten times the data statistics and is generated with a range of 
underlying T(D* + ) from to 130 keV. We do this for the three samples and compute offsets 
between the generated and fit values for the width and mean energy release. We also note 
that in all three simulated samples there are no trends in the difference between measured 
and generated width as a function of the generated width; the offset is consistent with zero 
as a function of the generated width of the D* + . Table [TJ summarizes this simulation study. 
We apply these offsets to the fit value that we obtain from the data. For the energy release 
all samples show small shifts, —7 ± 3 keV for the nominal, —12 ± 10 keV for the tracking, 
and —12 ± 5 keV for the kinematic. 

Figure ^ displays the fits to the three data sample. The results of the fits are summarized 
in Table |. Correlations among the floating parameters of the fit are negligible. 

The agreement is excellent among the three fits, and when the offsets from Table || are 
applied we obtain the results given in the last row. The uncertainties are only statistical. 

We discuss the sources of systematic uncertainties on our measurements of the width 
of the D* + in the order of their size. The most important contribution is the variation of 
the results as a function of the kinematic parameters of the D* + decay. The next most 
important contribution comes from any mismodeling of <tq's dependence on the kinematic 
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TABLE II. Summary of our data sample, simulation biases, and fit results. 





Sample 


Parameter 


Nominal 


Tracking 


Kinematic 


Candidates 
Background Fraction (%) 

Tfit — r genera ted (keV) 

Fit T (keV) 
D* + Width (keV) 


11496 
2.51 ±0.27 
2.7 ±2.1 
98.9 ± 4.0 
96.2 ±4.0 


368 
4.1 ± 1.9 
1.7±6.4 
106.0 ± 19.6 
104 ± 20 


3284 
4.05 ±0.49 
4.3 ±3.1 
108.1 ±5.9 
103.8 ±5.9 



parameters. We take into account correlations among the less well measured parameters of 
the fit, such as f m i. s and <J m is, by fixing each parameter at ±lcr from their central fit values, 
repeating the fit, and adding in quadrature the variation in the width of the D* + and Qo 
from their central values. We have studied in the simulation the sources of mismeasurement 
that give rise to smearing on the width of the D* + by replacing the measured values with 
the generated values for various kinematic parameters of the decay products. We have then 
compared these uncertainties with analytic expressions for the uncertainties. The only source 
of smearing that we cannot account for in this way is a small distortion of the kinematics of 
the event caused by the algorithm used to reconstruct the D° origin point described above. 
We have also checked that our simulation accurately models the line shape of other narrow 
resonances visible in our data. Notably the decay A — > pn~ , has a Q only seven times 
that of D* + — > D°tt^ ow . In the A decay we select the 7r~ to have a momentum in the 
range of those in the D* + decay, and the visible A widths agree to a few percent between 
data and simulation. We consider uncertainties from the background shape by allowing the 
coefficients of the background polynomial to float. Minor sources of uncertainty are from 
the width offsets derived from our simulation and given in Table [II], and our digitized data 
storage format. 

An extra and dominant source of uncertainty on Qo is the energy scale of our mea- 
surements. We evaluate this uncertainty by studying K s — > 7r + 7r~ decays in our data. In 
order to bring the K s mass central value in agreement with the nominal one, we make small 
relative momentum corrections, less than 0.3%, for tracks with momenta between 100 and 
500 MeV/c. These corrections only affect the slow pion. Applying these corrections to the 
momentum of the slow pion in our data we find a shift in the fit value of Qo, —4 keV for all 
the samples, and a negligible change in the width. We evaluate uncertainties in the energy 
scale by varying an overall momentum scale to change the K s — > 7r + 7r~ mass by ±30 keV, 
the uncertainty on that mass ||, and applying the statistical errors we obtain from the 
calculations of the momentum corrections discussed above. 

Table [TTT] summarizes the systematic uncertainties on the width of the D* + and Qo- 
In summary we have measured the width of the D* + by studying the distribution of the 
energy release in D* + — > D°n + followed by D° — > K~n + decay. With our estimate of the 
systematic uncertainties for each of the three samples being essentially the same we chose 
to report the result for the sample with the smallest statistical uncertainty, the minimally 
selected sample, and obtain 
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TABLE III. Systematic uncertainties on the width of the D* + and Qq 





Uncertainties in keV 
Sample 




Nominal 


Tracking 


Kinematic 


Source 


5T(D*+) 


5Qo 


5T(D*+) 


5Qo 


5T(D*+) 


SQo 


Dependence on Kinematics 


16 


8 


16 


8 


16 


8 


Mismodeling of <tq 


11 


< 1 


9 


4 


7 


< 1 


Hit" ( .nvrplflf lnnQ 


8 


3 


9 


4 


9 


5 


Vertex Reconstruction 


4 


2 


4 


2 


4 


2 


Background Shape 


4 


< 1 


2 


< 1 


2 


< 1 


Offset Correction 


2 


3 


6 


10 


3 


5 


Data Digitization 


1 


1 


1 


1 


1 


1 


Energy Scale 


1 


8 


1 


8 


1 


8 


Quadratic Sum 


22 


12 


22 


16 


20 


14 



T(D* + ) = 96 ± 4 ± 22 keV, (3) 

where the first uncertainty is statistical and the second is systematic. This is the first 
measurement of the width of the D* + , and it corresponds to a strong coupling [|I|] 

g = 0.59 ±0.01 ±0.07. (4) 

This is consistent with theoretical predictions based on HQET and relativistic quark models, 
but higher than predictions based on QCD sum rules. We also measure the mean value for 
the energy release in D* + — > D°n + decay 

Qo = 5842 ± 2 ± 12 keV, (5) 

where the first error is statistical and second is systematic. Combining this with the mass 
of the charged pion, 139.570 MeV with an uncertainty less than 1 keV [§], we calculate 

m D * {2 oio)+ - m D o = 145.412 ± 0.002 ± 0.012 MeV. (6) 

This agrees with the value from the Particle Data Group, 145.436±0.016 MeV, from a global 
fit of all flavors of D*-D mass differences. 
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